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Scattering by a Two Dimensional Groove in a Ground Plane 

Kasra Barkeshli and John L. Volakis 

\ 

\ 

[' 

Abstract 

| 

Higher order boundary conditions involve derivatives of the fields beyond the first 
and were recently shown to be more effective than traditional first order conditions in 
, modeling dielectric coatings and layers. In this report an application of a third order 

generalized boundary condition to scattering by a filled rectangular groove is 
! presented. Deficiencies of such higher order boundary conditions are addressed and 

a correction is proposed for the present case. As part of the process of examining and 

I 

! improving the accuracy of the proposed generalized boundary conditions, an exact 

i 

solution is developed and a comparison is provided with a solution based on the 

i 

standard impedance boundary condition. 




I. 


INTRODUCTION 


Traditionally, the standard impedance boundary condition (SIBC) [1] has been 
employed to simulate thin material layers on perfectly conducting objects. As is well 
known, however, the SIBC provides limited accuracy and is particularly applicable to 
lossy and/or high contrast dielectrics. This is primarily because it cannot model the 
polarization current components that are normal to the dielectric layer. As a result, the 
SIBC has been found to be best suited for near normal incidence, unless the coating's 
material properties are such that limit penetration within the material. 

The SIBC is a first order condition in that its definition involves a single normal 
derivative of the component of the field normal to the modeled surface. Recently[2], 
however, a class of boundary conditions were proposed whose major characteristic is 
the inclusion of higher order derivatives (along the direction of the surface normal) of 
the normal field components. These were originally introduced by Karp and Karal [3] 
and Wienstein [4] to simulate surface wave effects, but have been found to be rather 
general in nature. In fact, they can be employed to simulate any material profile with a 
suitable choice of the (constant) derivative coefficients. Appropriately, they are 
referred to a generalized impedance boundary conditions (GIBC) and can be written 
either in terms of tangential or normal derivatives provided a duality condition is 
satisfied [2]. Unlike the SIBC they offer several degrees of freedom and allow an 
accurate prediction of the surface reflected fields at oblique incidences. This was 
demonstrated in [2] for the infinite planar surface formed by a uniform dielectric layer 
on a ground plane. It was found that the maximum coating or layer thickness that can 
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be simulated accurately with a given GIBC was analogous to the highest order 
derivative included in the condition. 

The GIBCs offer several advantages in both asymptotic and numerical 
analysis of electromagnetic problems. For example, in the case of asymptotic/high 
frequency analysis, they allow an accurate replacement of a coating on a layer with a 
sheet boundary condition amenable to a Wiener-Hopf analysis [5,6] or some other 
function theoretic approach [7]. In numerical analysis, the profile of a coating can be 
replaced by a simple boundary condition on the surface of the coating. This eliminates 
a need for introducing unknown polarization currents within the coating or material 
layer and thus leading to a more efficient solution. 

In this report we examine a numerical application of a third order GIBC for 
scattering by a material filled groove in a ground plane. Since the GIBCs were derived 
for a coating without terminations, of particular interest in this study is the examination 
of their accuracy near those terminations. It is, unfortunately, found that they must be 
supplemented by more accurate conditions in the vicinity of material discontinuities. A 
procedure is, therefore, introduced that combines the exact and GIBC formulations. 
Since the exact solution is required for comparison purposes and in developing the 
hybrid formulation, it is presented in the first part of the report. This is followed by a 
discussion on the limitations of a formulation based on the SIBC. The integral 
equation based on a third order GIBC is presented next. This is solved by the 
conjugate gradient FFT method having an O(N) memory requirement. In contrast, the 
exact integral equation is not amenable to such a solution and must be solved by a 
matrix inversion approach having an 0(N 2 ) memory requirement. The hybrid 
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formulation is presented last and results are given which show the validity of the 
procedure as well as its limitations. 


II. EXACT SOLUTION 


A. Formulation 

Consider the H z polarized wave 


H 



(X -■!!!.) cos 4> 0 + ya n«> 0 


( 1 ) 


incident on the two dimensional cavity as shown in Fig. 1 and we are interested in 
computing the sattered field. The exact formulation for a slit in a ground plane has 
been developed in [8] and a similar formulation can be followed in the case of the 
groove. Below we develop an exact integral equation based on this procedure. 

In accordance with the equivalence principle the cavity may be closed by a 
perfect conductor and the equivalent magnetic current 


may then be introduced on the cavity's top surface, where E x is the x component of the 
total electric field at the aperture. Referring to Fig. 2 and applying continuity of the 
tangential electric field across the aperture, we have 
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( 2 ) 


Exn = Exn’ 
z M z = - z M' z 
M’ 2 = -M z 

In addition, continuity of the total tangential magnetic field across the aperture gives 
h| + H,° (M 2 ) = H, b (W 2 ) = H, b (-M 2 ) = - H, b (M 2 ) ; y = 0* (3) 

in which H t b denotes the tangential magnetic field at the cavity’s surface in region b and 
likewise H t ° is the tangential field at y = 0 + . 

From (3) an integral equation can be obtained after substitution of the 
appropriate field expressions. On the aperture (y=0), these are 

jk(x-^)oos0 o 

H,=2e 2 (4) 

w 

H, a (M z ) = - jk„Y 0 J 2 M z (x’) G(x, x') dx' 

0 

w 

= j M 2 (x’( H o a (k 0 1 x - x'l ) dx' (5) 

0 

in which Y 0 = 1 /Zq is the free space admittance and ko is the free space wavenumber 
and a factor of two has been introduced due to image theory. 

To find an expression for H t b (M z ) we require the cavity Green’s function. To 
find it we proceed as follows: Set 
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M = z M' z 
Frr'PZ 


( 6 ) 

(7) 


where *F is the scalar wave potential satisfying the two dimensional wave equation in 
the source free region 


f 2 2 'N 

9 4 . 5 + k 2 

, 2 3 2 +R b 'F = 0 in region b except y = 0. 

^dx dy j 

subject to the appropriate boundary conditions. Also 

(8) 

I = - VxFx^-'F-y^-'F 
3y 1 dx 

(9) 

H = -jk b Y b F 

(10) 


in which Y b = Y 0 jp , 1^ = k Q ys b p b and (e b , i^) are the relative constitutive 

parameters of the material filling the cavity. The scalar wave potential 'F can be 
determined by recalling the boundary conditions for the tangential electric field on the 
walls of the cavity. These are 


E x (x, 0) = - M z 


(ii) 

E x (x, -t) = 0 

- w\L.r° 

(12) 

Ey (0, y) = 0 
Ey (w, y) = 0 

> => ^-'F 1 =0. 

dX x-Qw 

(13) 


Choosing, 
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^(x,y)= ^ A p cos [kp(y +t )] cos-^- ; kp = \/ k b 2 - 

p-o 


fpn'P 
"w 


satisfies (12) and (13). The constants A p are then determined to satisfy (11). That is 


fr 1 ,r'£ A [> k p s in [M y+t) ] c ° s ^2-i = 

y y-° Jo y-° 


M, 


(14) 


C0S ^ = - M 2- 

p-0 


Multiplying both sides by cos and integrating, we have 


X A p k p sin V J C0S ^T cos ^ dx = - f f 
p-o 


M z (x) cos dx 


(15) 


yielding 


w 


prcx -i 

l n J I , 9 l*ly COS ... OX 

p wkpSinkpt I z ' ’ w 


A = 


i|M 2 (x) 


(16) 


in which 


s>= 


1 

2 


P = 0 
p > 0 


Therefore, 


H >z> = -iS,Y b £ 


p-0 


wkp sin k p t 


W 


| M z (/) 


cos 


^P7tX A 


V W j 


dx' 


cos 


K (y +t ) 


cos 


r pTO( N 

v^rJ 


(17) 
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B. Moment Method Solution 


Substituting (4), (5) and (17) into (3) yields the integral equation 


k 0 Y 0 


w k Y 

J M z (x') (k 0 | x - x'|) dx’ - j £ 

n P-0 


w 


T> P™ 

■cos 


k p tan (kpt) W 


J M; (JO 


pnx' , 
cos — dx' 


= 2 e 


jk(x-^-)cos<|) 0 


(18) 


To discretize the above integral equation we expand the magnetuc current in terms of 
pulses basis function as 


M z (x) 


= £c„M„ 

n-1 


(x) ; M n (x) = i 


1 x n -|-<x<x n + |- 

0 , else 


(19) 


and by substitution back into (18), we obtain 


k o Y o 


V — 
^ 2 


n-1 


O , 11V . . j k b Y b 


X- — 

2 


0 ogx-rfOtK-^g^jj^f 


- 4 • 


P KX ’ ^ , 

COS ^ | C 0 S ~W~ dX 

n-1 


v.A 
* 2 


= 2e 


j k (x - — ) oos «t> 0 


(20) 
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Employing point matching, this can now be written as 


[ Y mn] [CJ = [VJ 


( 21 ) 


where 


Y — = 


k 0 Y 0 


X„+ — 

''n 2 


mn 


f (2) 

J h 0 (K, l> 


X m - x'l) dx* - 


jkY 


b b 


2w 




£ 

p«0 



P* x ml 

e p cos 

W J 


k p tan (kpt) 


Xp+ — 
^ 2 


I 


COS 


^P7CX A 


w 


x„-- 


dx’ 

( 22 ) 


and 


V m = e 2 (23) 

with x m = , m = 1 , 2, 3, ... The elements of the square admittance matrix can be 

further simplified to 


1 2 
Y = Y + Y 

mn mn T 1 mn 


With 


mn 


and 


r 

l 


k.Y n 

O 0 


j2 k A 


H o a (Ko|x m -x„|) 


for n = m 


for n * m 


(24) 


(25a) 
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Y 2 = 

1 mn 


Vb 

2w 


I 

p~0 



f P* x ml 

£ p COS 

W J 


kptan (k p t) 


P^ x n P 7 ^ 

cos ... sine . 
w 2w 


(25b) 


in which y^ 0.5772 is the Euler constant. 


C. TE Scattering from a Narrow Groove 

In this section we examine the exact formulation and compare it to that 
obtained via application of the standard impedance boundary condition for a narrow 
groove. By rewriting the integral equation (18) as 

k W W 

■f- 1 < x '> h." <Ki i x • *i ) <* * w-ktn \ m dx ' 

0 b ° 0 


*J 


Vb f 

w 2-t 


cos 


P7tX 


w 


f kptan kpt W 


j m 2 M 


COS 


pnx* , jk(x-|)cos$ 0 

Ol-k *— 


w 


dx' = 2e 


(26) 


and setting ri b = jZ b tan (k b t) we obtain 


k 0 Y 0 


W 


/ 


M, 


(X') H if (k„ 


w 


| x - x’l ) dx’ + 




J 


M. (x’) dx’ = 2e 


jk(x-^)<ns<!> 0 


(27) 


0 0 

where we have included only the dominant mode (p=0) in the representation of the 
cavitiy’s Green's function. If we further assume a constant tangential electric field 
variation (M z (x) ~ constant) over the aperture, (27) reduces to 


k 0 y 0 

2 



(2 1 jk(x-^)cos<D 0 

H 0 (k 0 I x - x’l ) dx' + ^- M z (x) = 2e 2 


(28) 
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i which is the integral equation resulting from an application of the standard impedance 

i 

boundary condition (SIBC). 

‘ The approximations required to go from (26) to (28) are not obviously 

expected to hold unless the groove is small in depth and width, say k<,w < 0.6 and 

i 

j k<jt < 0.6. Consequently one does not expect the SIBC to provide a good simulation of 

the cavity scattering, particularly near grazing incidences where M 2 (x) in not symmetric 
and may have rapid spatial variation. In fact, it may also be observed that the SIBC 
integral equation , at best, yields the average of the actual equivalent current 
distribution on the surface of the crack. This is clearly illustrated in Fig. 3 where the 
currents and echowidth predicted via (28) and the exact integral equation (20) are 
compared. The result shown in Fig. 3 is a typical situation and it is not surprising that 
the backscatter echowidth predicted by the SIBC formulation is in substantial error at 
oblique incidences. 

III. FORMULATION WITH GENERALIZED IMPEDANCE BOUNDARY 
CONDITIONS 

It is desirable to work with a formulation (or an integral equation) that is 
amenable to a conjugate gradient FFT (CGFFT) implementation. The CGFFT has an 
0(Nf memory requirement and can thus be suitable for treating large size grooves or 
cavities, particularly when applied to three dimensional geometries. Unfortunately, the 
! exact formulation, in addition to being restricted to rectangular grooves and cracks, is 
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not suitable for a CGFFT implementation. On the other hand the integral equation (20) 
resulting from an application of the SIBC, a first order condition, although amenable to 
a CGFFT is not of acceptable accuracy. Recently, however, higher order impedance 
boundary conditions involving field derivatives beyond the first have been found to 
provide a substantially better simulation for fairly thick dielectric coatings. These are 
referred to as generalized impedance boundary conditions (GIBCs) and take the form 



(30a) 


1 


a ’m 


a (-jk 0 ) m ay m 


= 0 


(30b) 


where a m and a' m are constants specific to the surface, layer or coating being 
modeled. For M=1 , they reduce to the SIBC provided we set 


Wo 



tan (ly) = 




( 31 ) 


where N = . A third order GIBC corresponds to M=3. In that case, an accurate 

simulation of the reflection coefficient for a metal-backed uniform dielectric layer can 
be obtained by choosing 


1 1 


a o-( N ' 2N^ 


tan (ktN) - tan ) 


a i = le b 


1 + tan (ktN) tan (^-) 


( 32 ) 


a 2 = 2 n I tan (ktN) - tan + kt (N - ^) 


1 + tan (ktN) tan (^) 


ikte, 


a 3 ' 2N 


tan (ktN) - tan (^ ) 


and 


a' =(2N -1) 


1 + cos (ktN) cot (-^j) 


a ’i = 1 2Np b 


cot (ktN) - cot (^) 


(33) 


a ’ 2 = 1 + cot (ktN) cot (-^) + kt (N - ^j-) 


cot (ktN) - cot (^j-) 


--"^L 


1 + cot (ktN) cot (J^) 


The above conditions are applied on the surface of the coating and predict the proper 
surface wave modes. However, they were derived for an infinite layer without the 
presence of any terminations. Therefore, when applied to the case of a groove having 
abrupt material terminations at x=0 and x=w, we expect that the simulation provided by 
(30) in conjunction with (32) and (33) will not be as accurate. As a result, the GIBC 
must be supplemented by additional (more accurate) conditions at the terminations of 
the coating or in this case the groove. At this point, no standard methodology has 

1 2 



been devised for imposing these supplementary conditions, but in either case such 
conditions will be specific to the geometrical and material properties of the termination. 
Before, however, we examine the issue of supplementary conditions, let us first 
proceed with a direct implementation of the third order GIBC noted above. It will be 
seen that comparison of the results obtained via a direct application of this GIBC will 
guide us on how these can be supplemented at the terminations. 

For H z polarization, H y =0, and thus the relevant GIBC is (30a). Expanding this 

we have 


( 


7f 

3 n*"T -2-r- 

k dy 


jk 


jk° 3y 


3y 


= 0 ; y = 0 


(34) 


To introduce the equivalent magnetic current M z = E x in (34), we note that 


V- E = 0 


(35) 


and thus 


dE y 3E X 3M Z 
By ~ " Bx = " Bx 


(36) 


Substituting (36) into (34) now yields 


a 2 

a + ~ 

0 k 2 
K o 


k 2 + -^- 
0 Bx 2 


W 


1 


(x') ^ Hf (k | x - x'| ) M * 


J 0 


13 


f 


J k c 


a + — : - 
1 k 2 

K o 


.2 

k ° + TT 
dx 


3M Z (x) 
9x 


f 


a 


V 


K ay 


\ 




2Z n cos H z (37) 


In deriving (37) we also employed the wave equation 


% 3 % 2 

ay 2 " ax 2 ' 0 


and have set 


(38) 


Ey-Ejj + Ey' + Ey 8 - - 2 Z,, COS ^ H* + 


(39) 


where 


W 

s i r ( 2 ) 

E y -ij2M 2 (x') H 0 (kj, | x - x'| ) dx' (40) 

0 

is the y component of the scattered field in which the factor of 2 is due to image theory. 

Integrating both sides of (37) with respect to x eliminates one of the 
derivatives. Doing so, we obtain 

W 

k, Y„ e (2i 

J M z (x')H q (k 0 |x - x’| ) dx’ + 

o 


0 


ni-L 

K, 3x 


+ 


( 



1 + 




X M z (x) = 2 







j k( x - — ) cos 

e 2 


J 


♦o 


( 41 ) 


and we note that for a 2 = a 3 = 0 (41 ) reduces to the SIBC integral equation (28) 
provided the identification noted in (31) is also employed. 

The integral equation (41) derived by imposing the GIBC lends itself to a 
solution via the CGFFT method. Defining the fourier transform pair 

oo 

G (k x ) = / { G ( x )} = J G(x) e JlV< dx (42) 

- oo 
oo 

G(x) = / ' {<3 og} = | G (X) e' M dk, (43) 

“OO 

we have 


/ h" < k o M>} = H ? (k 0 |X|) = -jA=f = — % 


VkTC k/f x 2 -1 


(44) 


(/-T = j) and by recalling the convolution theorem 


W 

f M 2 (x') H 0 ' 3 (K, I X - Xl ) dx' = 2j JhM. 

0 K.V‘x-1 


(45) 


since M z (x) is zero outside 0 < x < w. Also, in the transform domain 



_ 1 _ 

Kf 3x2 



and (41 ) may thus be written as 
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(46) 



A 


M Z (U 


I jk(x-^)cos<|> 

f = 2Z»e 2 


0 


A CGFFT implementation of (46) is a straightforward task and typical results as 
computed by (46) are shown in fig. 4. As expected, the third order GIBC applied to the 
groove, predicts reasonably well the exact magnitude and phase of the current 
distribution when away from the termination of the groove. The accuracy of this 
prediction, of course, depends on the depth of the groove, t, and the material 
properties of the dielectric filling. Our preliminary investigation indicates that for 
lossless dielectrics, the above third order GIBC provides a reasonable prediction of the 
current distribution away from the groove terminations for kt < 1 . However, for lossy 
dielectrics substantially deeper grooves can be simulated. 

Next, we consider a hybrid GIBC-exact formulation to alleviate the difficulties 
of the GIBC in predicting the currents near the groove terminations. 


IV. HYBRID GIBC-EXACT FORMULATION 


The GIBC formulation in conjunction with the CGFFT solution method offers 
the substantial advantage of having an O(N) menory requirement. However, as seen 
in Fig. 4, the current distribution predicted by the third order GIBC is not of acceptable 
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accuracy when within 0.2 wavelengths of the groove terminations or so. To alleviate 
this difficulty, one approach is to feed the currents predicted by the GIBC integral 
equation (41) away from the edges into the exact integral equation (20). The last can 
then be solved for the remaining currents in the vicinity of the groove terminations. 
This only requires the inversion of a small matrix thus alleviating the usual difficulties 
with storage. 

Q 

Suppose now that M z (x) denotes the current computed via the GIBC integral 

equation given in (41) and likewise M®(x) denotes the current computed via the exact 

G 

integral equation. Employing M z (x) in place of M z (x) in (20) for x A < x < w - x A yields 


XqYq , 


"a w 

J M ® (x') (k 0 |x - x’l) dx' + J M z ® (x*) (k Q |x - x’|) dx’ 


W-X A 


-J 


■ ki>Yb 2 


r 


P ,nc P7tX 

cos 


w ^V an (V) w ^ 


<{ jM®(x’)cos^-dx’ + 


w , 1 

Jm®(x') cos dx' ^ 

J 


W-X A 


= 2 e 


jk(x-^)oos<|> 0 ■ ” A 


k o Y o f t» Q u (2) 
2 


J M z (X 1 ) Ho (k 0 1 x - x 'l ) dx ’ 




W- X A 


2 

p-0 


P” x f - 
w Lu kptan(kpt) W J z [ ' 


cos 


P7CX' . , 
cos dx’ 


(47) 
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Q 

Assuming that M z (x) has already been determined via a CGFFT solution of (46), the 
entire right hand side of (47) is known. Thus for x A £ 0.25 a 4x4 or a 6x6 square 
impedance matrix is required for the solution of M z (x). In general, continuity of the 

Q 6 

current density must also be imposed at the transition regions between M z (x) and M z (x). 

The results shown in figure 5 clearly show that the proposed hybrid 
formulation can provide an accurate prediction of the scattering by a groove. The 
bistatic and backscatter echowidths presented in the course of the above 
developments have been summarized in Fig. 6. Other examples for the current 
distributions and corresponding echowidths are given in figures 7-11. In these 
figures the following labeling has been employed 

EXACT : Data from a numerical implementation of the exact integral 

equation (20). 

SIBC: Data from a CGFFT implementation of (46) with a$ = a 2 = 0 and in 

conjunction with (31) 
condition. 

GIBC-3: Data from a CGFFT implementation of the integral equation (46) 

resulting from the third order generalized impedance boundary 
condition. 

Hybrid-1 : Data from the hybrid SIBC-exact formulation. 

Hybrid-3: Data from the hybrid GIBC-exact formulation employing the 3rd 

order GIBC. 
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V. CONCLUDING REMARKS 


An application of a third order generalized boundary condition (GIBC) to 
scattering by a two dimensional dielectrically filled cavity was considered. In the 
process of examining the accuracy of the GIBC, an exact solution was developed and 
a solution based on the standard impedance boundary condition (SIBC) was 
examined. An analytical comparison of the integral equation based on the SIBC with 
the exact, revealed the well known limitations of the SIBC formulation. It was 
concluded that the SIBC integral equation will, at most, generate an average of the 
actual current distribution provided the groove is very shallow. 

The GIBC integral equation was found easier to implement. Furthermore, 
unlike the exact integral equation, it was amenable to a conjugate gradient FFT 
solution and is, thus, attractive for three dimensional implementations. It was found to 
predict the correct current behavior reasonably well away from the terminations of the 
groove particularly for lossy dielectric fillings. However, the inadequacy of the GIBC 
formulation near the groove terminations proved problematic. The GIBC conditions 
needed supplementation in these regions and several approaches were examined to 
correct their deficiency. Our initial hope was that the addition of filamentary currents at 
the edges would provide the required correction as was already done in the case of an 
isolated thin dielectric layer. This approach, however, was not found suitable for the 
subject geometry. Instead, the incorrect currents near the groove terminations were 
replaced with those computed via the exact integral equations. Specifically, the 
currents computed via the GIBC formulation away from the groove termination were 
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employed in the exact integral equation to generate a small 4x4 or a 6x6 matrix for the 
currents in the vicinity of the terminations. This was referred to as a hybrid exact-GIBC 
approach and was found to provide a reasonably good simulation of lossy dielectric 
fillings at all angles of incidence and observation. In case of lossless and low contrast 
dielectrics, the simulation was adequate for groove depths up to 3/20 of a wavelength. 
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Fig. 1 . Geometry of the rectangular groove in a ground plane 
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Fig. 2. Illustration of the application of equivalence and 
image theory. 
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Fig. 3 Comparison of results for scattering from a groove obtained by the 
exact and the SIBC formulations; groove width w-lA, depth WJ.2A, e^.O-jl.0, 1. 






Fig. 4 Comparison of results for scattering from a groove obtained by the 
exact and the third order GIBC formulations; groove width w=lA, depth 
fc=0.2A, e k =7.0-jl.0, /z k = 1. 


i 







Bistatic Echo width/A . dB | Magnetic Current Density |M,(x)/E,| 


Groove Current Distribution 


Groove Current Distribution 



0.00 0.25 0.50 0.75 1.00 


Aperture Position. x/X 




Backscattering Pattern 



90 120 150 180 


Angle of Incidence. p,deg. 


Pig. 5 Comparison of results for scattering from a groove obtained by the 
exact and the hybrid GIBC formulations; groove width w-lA, depthtwOJZA. e k -7.0-jl.0, #* k -l. 
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Fig. 6(a) Comparison of bistatic scattering patterns obtained by different 
formulations as indicated; w=lA, b=0.2A, e h =?7.0-jl.0, /u b =l. 
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Fig. 6(b) Comparison of backscatterinc patterns obtained by different 
formulations as indicated; w=lX, tr=0.2X, c b =7.0-jl.0, /x h =l. 
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Fig. 7 Comparison of results for scattering from a groove obtained by 
the exact, third order GIBC alone and hybrid formulations; groove width 
w=lA, depth tF0.05A, c b =l, /z b =l. 
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Fig. 8 Comparison of results for scattering from a groove obtained by 
the exact, third order GIBC alone and hybrid formulations; groove width 
w=lA, depth fc=0.4A, c fc =^7.0-jl.0, /z h =l. 
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Fig. 9 Comparison of backscattering patterns obtained by different 
formulations as indicated; w=lA, t=0.4A. c h =7.0-jl.0, /Lt b =l. 
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